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Abstract 

We consider a pair of parallel D p _i and anti-D p _i branes in flat space, 
with a finite separation d along some perpendicular spatial direction and at 
finite temperature. If this spatial direction is compactified on a circle then 
by T-duality, the system is equivalent to a D p -antiD p pair wrapped around 
the dual circle with a constant Wilson line A « d on one of the branes. We 
focus in particular on the p = 9 case and compute the free energy of this 
system and study the occurrence of second order phase transitions as both 
the temperature and Wilson line (brane-antibrane separation) are varied. In 
the limit of vanishing Wilson line we recover the previous results obtained in 
the literature, whereby the open string vacuum at the origin of the tachyon 
field T = is stabilized at sufficiently high temperature at which a second 
order phase transition occurs. For sufficiently large Wilson line, we find new 
second order phase transitions corresponding to the existence of two minima in 
the tachyon effective potential at finite temperature and tachyon field value. 
Entropic arguments suggest that as the system cools, the tachyon is likely 
to find itself in the minimum that approaches infinity as the temperature 
vanishes (i.e. the one corresponding to the closed string vacuum), rather than 
the minimum at T = (corresponding to the open string vacuum). 
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1 Introduction 

The study of unstable (non-BPS) D-brane configurations in flat space [1] has been a 
fertile area of research in recent years. Sen's conjectures |T] concerning what happens to 
unstable D-branes and the fate of the open string vacuum has been supported by results 
in boundary string field theory, (BSFT) [2], [3], pE]. In this picture a coincident parallel 
brane-antibrane configuration is unstable to decay through the open string tachyon field 
T rolling down to its minimum and thus producing enough negative energy to cancel that 
coming from the brane tensions. Thus the open string vacuum decays to that of the closed 
string. A comprehensive review of tachyon condensation from the point of view of BSFT 
and other approaches can be found in [5]. 

Whilst the above behaviour of the open string tachyon is true for a system at zero 
temperature there have been several papers discussing the situation if the brane-antibrane 
pair is considered as part of a thermodynamic system at finite temperature [6], [7j, [8], [9], 
|10j . Including finite temperature effects is interesting because it's possible that such brane 
configurations could survive in the early universe and thus be stable at finite temperature 



In [7j, Hotta investigated the phase structure of a finite temperature D p brane- ant i D p 
brane pair (which we will abbreviate in this paper as D p — D p ) where the branes were 
assumed to be coincident and in flat space. Using the framework of boundary string field 
theory (BSFT) he showed that in the p = 9 phase transition occurs just below 

the string Hagedorn temperature, whereas for p < 9 there is no phase transition. In the 



en- 
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p = 9 case, the zero temperature minimum of the tachyon effective potential was shown 
to shift from T — > oo towards T — > as the temperature approached criticality. Thus the 
interpretation is that the open string vacuum is stabilized at sufficiently high temperature 
(but below the Hagedorn transition) in the case of D 9 — D 9 whereas for pairs of lower 
dimensionality, no such transition occurs and the point T = remains unstable at high 
temperature. 

These results are in broad agreement with those of Danielsson et al [6] who investigated 
the same system but rather than including the full set of string states, they focussed on 
the truncation to the tachyonic sector only, in computing the free energy. 

In this paper we wish to generalize the results above to the case where the D p _i — D p _i 
pair is separated along some perpendicular spatial direction, but still parallel and in flat 
space. We shall assume that the pair has a finite separation d along a perpendicular 
spatial direction which is compactified on a circle S . Then by T-duality, the system is 
equivalent to a D p — D p pair wrapped around the dual circle S 1 , with a constant Wilson 
line A « d turned on one of the branes [9J , [H] . 

At zero temperature, one may extend the BSFT results to include Wilson lines and 
obtain an expression for the effective potential at 1-loop V e ff(T, A) depending on T and 
the Wilson line A. At tree level, the extrema of this potential depend on the size of A. 
For A < -7= the potential has a local maximum at T = as in the case of a coincident 

V2ct' 

brane-antibrane. If A > T = becomes a local minimum and so the open string 

vacuum is metastable [5J. We shall see that when we consider this latter situation at 
finite temperature, we have an interesting situation whereby the effective potential (in the 
canonical ensemble) has two local minima at finite values of T, which then approach the 
values T = 0, oo, as the temperature approaches zero. Thus we can ask the question what 
vacuum the tachyon field will likely be found, i.e., which vacuum is thermodynamically 
favoured over the other? Also what is the likelihood of a first order phase transition 
occurring since we anticipate that the two minima might well become degenerate at some 
particular temperature, so that quantum tunnelling may become important. 

The structure of this paper is as follows. In section 2 we review the two derivative 
truncation of the BSFT approach to studying the tachyon potential for coincident D p _i — 
Dp-i system [2], [3] and its extension to the case of finite separation (or addition of 
Wilson lines being included in the world-volume action of the T-dual wrapped D p — D p 
system). In section 3 we consider the 1-loop (annulus or cylinder) computations of the 
free energy of open strings stretched between separated D p _i — D p _i pair. In section 4 
we then investigate the critical points of the free energy and determine the nature of the 
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phase transitions as the temperature approaches the Hagedorn temperature from below. 
In particular we compare the situation when the Wilson line modulus A is greater than 
or less than its critical value A cr u = -7= . Finally in section 5 we draw some conclusions 
from our results. 



2 Two- derivative truncation of the BSFT of the DD 

system 



In string theory a pair of parallel D p — D p pair constitutes an unstable object. 
To study the dynamics of unstable D-branes, the BSFT [2], [3] is a useful tool and it 
has provided a good understanding of tachyon condensation at the classical level. It 
describes the off-shell dynamics of open strings in a fixed on-shell background of closed 
strings in which an open string field configuration corresponds to a boundary term in the 
world-sheet action of the string. Therefore, specifying a boundary term means giving the 
background values of the various modes of the open string. It is based on the Batalin- 
Vilkovisky formalism whose master equation provides the effective action of the theory. 
In the bosonic string theory, the disk partition function of the open string theory Z and 
the BSFT action are related by the master equation 



where g l are the couplings of the boundary interactions and (3 l are the corresponding 
world-sheet /3-functions. Given a specific form of the tachyon profile, the BSFT action 
reduces to the effective action for the tachyon field allowing us to compute the tree level 
tachyon potential. For superstrings, the tachyon /^-function is zero and eq. (12.11) reduces 
to 



The partition function Z was computed in [T3]. Let us briefly review their results here. 
The disc partition function is formally defined as 




(2.1) 



S = Z 



(2.2) 




(2.3) 



where 




(2.4) 



is the bulk action for the NSR string. 
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The boundary term of the D p — D p system is computed introducing auxiliary boundary 



fermion superfields T 1 



OF 1 where I = 1, 2 m , and, N 



im-1 



is the number of 



Dp — Dp pairs. Consider, for example, the case where we have 2 m branes. The 2 m x 2 m 
matrices of the gauge group U(2 m ), generated by the branes, can be expanded in terms 
of SO(2m) gamma matrices. Now, instead of gamma matrices, one can introduce 2m 
boundary fermion superfields Y 1 with action S = — f drdO^T 1 DT 1 , and, after canonically 
quantizing, one arrives at the anti-commutation relations {t] 1 ,^} = 2S IJ . Thus, t] 1 can 
represent the Clifford algebra needed for the expansion of the 2 m x 2 m matrices [T3] . 
In the case, e.g., of a single Dg — Dg pair, expanding the resulting action in terms of the 
component fields one has 



Si 



bndy 



1 

4' 



a 
~2 



dr 



(2.5) 



Here I,J = 1,2, 



.4 



± 

D..T 1 



d^T 1 - iA^T 3 



(2.6) 



and the gauge fields A^ on the brane and anti-brane, respectively, have been expressed 
in terms of the abelian gauge fields Aj/ , (anti-symmetrized in J, J) and A^. Moreover, 
Fftv = d [lx A v ] and F IJ p V = d^A 1 - 1 ^. 

In the case of a constant tachyon field and zero gauge fields, the boundary action 
reduces to 



5, 



a 



bdry 



drTT 



Irpl 



(2.7) 



Since there is no other dependence on the tachyon field in the bulk action, we learn that 
in this case the tachyon potential for the Dg — Dg system is 



Vo0n = 2T 9 e- w l T l 2 



{21 



where we defined T — ~ (T 1 + iT 2 ), whereas Tg denotes the tension of a Dg-brane which 
is defined, for general p, by 

T P = - - \ P+1 (2.9) 
(2ir)Pa'— g s 

where g s is the string coupling constant. The stable vacuum is at T = oo, where the 
vacuum energy vanishes. Since the potential (12.81) is exact, it gives a proof of Sen's 
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conjecture [T] that the negative energy contribution from the tachyon precisely cancels 
the D-brane tension: under tachyon condensation, the D-brane will decay into the closed 
string vacuum without any D-branes, therefore, excitations are described by closed strings 
alone. 

Let us turn now to the case of a spatially dependent tachyon. In this case, by a 
combination of spacetime and gauge rotations one can bring T 1 to the following form: 

Va'T 1 = u 7 X 7 (2.10) 

where u 1 are constants. When the gauge fields are zero, one can compute the partition 
function Q using eqs. (J23D, and (12TTUD . The result is [T3] 

r - 2 

Z = 2T 9 / ftoe-^fjFfW^/) 2 ) (2.11) 

J i=i 

where 

= 4 XL { X) 

V ; 2r(2x) V ' 

The partition function allows us to have an expression for the action of the tachyon at all 
orders in derivatives. However, there is an ambiguity in the expansion, because any term 
with at least two derivatives acting on T can be added. At quadratic order the result (for 
the case of coincident D 9 — D 9 ) is unambiguous 

S^2T 9 J d l0 x e~ WTf [1 + 8W 2 ln(2) d"T d^T + . . .] (2. 13) 

where the expansion 

F(x) = 1 + 21n(2)x + C(x 2 ), x -> (2.14) 

has been used. 

Now, consider the case where one of the spatial directions, y, is wrapped on a circle 
of radius R < \fa' and that we have a constant Wilson line A wrapping the compact 
direction on say the D 9 brane. The gauge field strength in (12. 5p vanishes and the only 
dependence on the gauge field comes from the covariant derivative. We can lift the above 
expression to include the covariant derivative by simply changing the argument of the 
function F. 

Applying a T-duality transformation along y, the gauge field is mapped to the Higgs 



1 In a recent paper |14j . the authors have investigated tachyon condensation in the separated D p — D p 
system including the effects of higher level terms, in supcrstring field theory. 
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A<A, 



A>A C1 




Figure 1: Left: Tachyon potential for < A < A cr . Right: A > . In all the plots 
T 9 = 0.1, a' = 1/2. 



field which measures the distance d between a D 8 — D 8 pair, separated along the dual 
coordinate y with d ~ \A\. 

Adopting the normalization of the tachyon field used in [5] , the action (12.131) becomes 



S = 2T 9 J d 9 xdye~\ T \ 2 [l + 2a'\d ll T\ 2 + 2a'A 2 \T\ 2 ] 



The potential term is 

The extrema are given by 
dV (T) 

d\T\ 



V (T) = 2T 9 e~^ 2 [l + 2oM 2 |T| 2 ] 



2T 9 |T|e-l T l 2 (4A 2 a' - 2 (2A 2 a'\T\ 2 + l)) = 



i.e., 



\T\ = 0, \T\ = +oo, and \T\ 



V2A* 



a' 



[2-15) 



(2.16) 



(2.17) 



(2.18) 



V2a'A 

To study the nature of these extrema, we need to compute the second derivative: around 
\T\ = we have 

d 2 V (T) ] 



d 2 \T\ 



| m=0 = m 2 = 4T 9 (2a' A 2 -1) 



(2.19) 



Therefore, we see that this potential has a minimum at \T\ = if A > y or it has a 

true tachyonic instability if A < \J^r- Figure [T] shows the different cases. 
This behavior has a clear physical interpretation: recall that our model is equivalent to the 
case of a parallel D§ — D§ pair separated by a distance d. If the distance d is large enough, 
then the tachyon mode between the two should go away, since the tachyon field comes 
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from the open string suspended between the two branes and thus that string acquires a 
mass lift when two branes are distant. 

Notice also that in order to get a canonical kinetic term in the BSFT action we must 
perform the following redefinition of the tachyon field: T = T(0) with 

r\T\ 

= v^a'Tg / dse~ s2/2 (2.20) 



Jo 

With this redefinition, the action (12.151) becomes 

S = J d\dy Q(90) 2 + K>(T(0))) (2.21) 

and the tachyon vacuum at infinity is placed at a finite value of the new field 0. Indeed, 
the two local minima are 



'0 



0, fa = y/Ana'Tg. (2.22) 



This redefinition allows us to compute the mass of the tachyon: in the presence of a 
Wilson line A it is given by 



M2 = an M = 1 

dd) 2 a' 



Tf - 1 
2 



+ a'A 2 (|T| 4 -4|T| 2 + l) 



(2.23) 



whereas if A = we have 

1 ITI 2 — 1 

ML„ = iX_ (2.24) 



Notice that the same results were found in [15] but with different methods. 
Henceforth, we will consider only the real part of the tachyon field: this is consistent with 
the tachyon equations of motion and it is also a natural setup since we are not interested 
in lower dimensional D-brane left after the tachyon condensation which needs complex 
tachyon configurations. 

3 Free energy of open strings stretching between a 
Dp — D p pair 

Before we discuss the free energy of strings stretched between a D p — D p pair, let us first 
comment on the issue, raised by Hotta in j7], concerning the microcanonical ensemble 
vs canonical ensemble framework for the computation of the tachyon finite temperature 
potential. It was shown in [7] that whilst in principle, the microcanonical picture is more 
trustworthy as we approach the Hagedorn temperature, in fact for the case of a coincident 
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Dg — Dg pair, the micro and canonical ensembles agree in the nature and existence of the 
second order phase transition of the tachyon effective potential near the origin. For the 
case p < 9 the case is less clear as the predictions for phase transitions do not entirely 
overlap for the various values of p in the two formalisms. In this case it is better to adopt 
the microcanonical ensemble as in [TJ, [8]. 

Since we wish to consider the case where we turn on constant Wilson lines around 
a compact spatial S l in the D p — D p system, we should consider how this affects the 
predictions made in both frameworks. In fact in all cases, the additional terms in the free 
energy (the singular part of which is used in [7J, [8] to extract the density of states in the 
microcanonical ensemble) coming from the Wilson line A can be computed. As we shall 
see below, the Wilson line only appears as an effective shift in the tachyon mass term 
when one considers the sums over all states contributing to the 1-loop partition. As such 
one can verify that at least in the p = 9 case, the canonical and microcanonical formalism 
will agree as regards the nature of the phase transitions in the tachyon effective potential 
even with A ^ 0. For p < 9 and i ^ one should again adopt the microcanonical 
ensemble too. It is straightforward to extend the techniques and results in [7j to include 
A but for brevity we will simply use the canonical ensemble in this paper, in which the 
one-loop part of the tachyon effective potential is given by the free energy of open strings. 

Thus we will primarily focus on the p = 9 case in this paper and do all computations 
in the canonical ensemble. By T-duality this is equivalent to separated D 8 — Dg pair with 
separation d ~ A. It is interesting to see whether the finite temperature could drastically 
modify the tachyon potential: in the case of zero Wilson line this is motivated by the fact 
that the tachyon field at T = can become stable and there is no tachyon condensation 
[BJ, [7J. In the presence of a brane separation, we might also be interested in the fate of 
the metastable minimum at T = (see Figured]). 

Temperature corrections to the potential (I2.16P come from the evaluation of the path 
integral ( 12. 31) over all connected graphs of strings on the space where the Euclidean time 
direction is compactified on a circle of radius equal to the inverse of the temperature (3. 
We will consider the weak coupling approximation in which the strings can be thought as 
an ideal gas, that is to say, ignoring the interactions of open strings. We take into account 
only one-loop amplitude considering only zero-genus oriented Riemann surfaces. 

Let us first quickly review the situation for a coincident D p — D p pair (though following 
our discussion above we will ultimately focus on the case p — 9). The effective potential 
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at finite temperature is given by 



V eff (T,P) = V (T) + V 1 (T,P) 



(3.1) 



where Vi(T, j3) is the one-loop finite temperature potential. 

Since we work in the canonical ensemble the one-loop part of the effective potential above 
is related to the free energy F(T, (3) of open strings: 



where V is the volume of the system, the D p — D p pair in our case. 

At this point, one immediately faces difficulties: in order to compute V\_(T,j3) or F, 
we need to include quantum corrections to the BSFT. At tree level the BSFT action 
is essentially given by the partition function on the disk and at one loop one might 
expect that the first loop correction corresponds to the partition function on a world- 
sheet of cylinder or annulus topology. However, because the boundary interactions break 
conformal invariance this result would depend on the choice of the Weyl factor. 
Nevertheless, there have been several attempts to generalize the BSFT to the one-loop 
amplitude in the D p — D p system [19] . All of them assume that the relation ( 12.2ft is still 
true at one loop. Then, they construct the partition function at one loop by keeping fixed 
the boundary of the disk and the tachyon profile on it and adding more boundaries and 
handles to the string world-sheet diagram. In particular, on the annulus one has 



where u is the coefficient of the linear tachyon profile eq. (12.101) . a nxc d is the boundary of 
the disk and b is the inner boundary of the annulus. Similarly, the cylinder amplitude can 
be computed in the closed string channel using the boundary state formalism. It is well 
known that the partition functions obtained in the two different schemes agree on-shell 
thanks to the open-closed string duality. However, the presence of the tachyon takes the 
theory off-shell and it is not clear, a priori, that the two different schemes yield the same 
result. In particular, since the boundary interactions are due to non- primary fields, the 
use of conformal maps to transform one worldsheet into another one is not helpful because 
the transformation laws of the fields are unknown. However, it seems that at one-loop at 
least, the two results are equivalent. 

In [16], for example, the partition function on the annulus and cylinder were computed 
in the presence of a constant tachyon profile. To fix the problems coming from the breaking 



V eff (T,[3) = V (T) + V- 1 F(T,P) 



(3.2) 




annulus 




(3.3) 
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of conformal invariance, they proposed to use a comparison with field theory results |17j . 
[18]. Indeed, given the partition function, one can in principle extract the contribution 
due to the tachyon and fix the background by comparison with the corresponding field 
theory results computed from the tree level effective action. This leads to equivalent 
expressions for the partition functions computed for the annulus and cylinder. 
The one-loop amplitude on the cylinder in such background is given by 



16 1 71% 



dr 



(Atxt) 



X 



9 3 (0\iT) 
6[(0\ir) 



IT 



6[(0\ir) 



(3.4) 



where V p is the volume of the _D p -brane. 

We would like to extend this result in order to include a more general background, i.e., 
a background in which the tachyon has a linear profile (I2.10p and to the case where the 
D p _i — D p _i pair are separated along a compact direction (equivalently turning on a 
Wilson line on the dual circle wrapped by a D p — D p pair) An easy way to do this comes 
from the following observation. 

In [7], it was noted that the one-loop amplitude (13.41) can be obtained by considering the 
free energy of strings stretching between the coincident D p — D p pair \ 



dr 



(Attt) 



oo 

P+l V — V \ ^ 

2 Z^2^ ex P 



M 2 r=l 



r 2 8 2 

-2na'M 2 NS r - tt-j- 



+ 



Yn 



(27ra')^ 



exp 



M R T=1 



-1t\olM 2 r t 



P 2 H t 

r 2 p 2 
fa 



7T 



(3.5) 



with the following mass spectrum 



M 2 NS 



M 2 R 



1 

~' 
1 



N B + N NS + 



T- 



N B + N R + 



2 

rp 2 

~2 



(3.6) 
(3.7) 



where M^s and Mr are the masses of the Neveu-Schwarz and Ramond sectors, respec- 
tively, whereas Nb, Nns and Nr are the oscillation modes of the bosons, Neveu-Schwarz 
fermion and Ramond fermions. Notice that the lowest mode of the NS sector (13.61) 
coincides with the mass of the tachyon field (I2.24p of the coincident D p — D p pair. 

2 We adopt the following definition for the Hagedorn temperature 0jj = 87r 2 a'. 
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This suggests a straightforward generalization of eqs. (13.61) and (13. 7p to the case of 
separated D p _i — D p _i. The only difference with the case described above is that in 
our model we have a constant Wilson line turned on on a circle of radius close to the 
string scale. Therefore, in general, we have to include quantized momenta in the direction 
parallel to the D p — D p system, and winding modes in the direction transverse to it. As 
for the presence of the Wilson line, notice that in the T-dual picture, the dependence on 
the constant Wilson line A in the tachyon mass (12.231) factorizes out, so we require that 
the lowest mode of the NS sector coincides with the tachyon mass. 
In the general case of D toroidal-compactified directions and d non-compact ones, the 
mass spectrum is given by [8] 



p—d 



M 2 NS 



i=i 

p—d 



M 



R 



E 



i=i 

where we have defined 



TTlj 

~R~i 



D 

+ E 

i=p— d+1 
D 

+ E 

i=p—d+l 



rijRi 



a' 



a 1 



rp 2 

2 



1 

2 



-\N B + N NS + - - - 

cr 

1 



Mi 



T' 



N B + N R + ^- +Ml 



Mi 



(3.8) 
(3.9) 

(3.10) 



- a A 2 (T 4 - 4T 2 + l) 

Inserting these two expressions into eq. (13. 5p and expressing the sums in terms of the 
^-functions using the conventions of [7j, the free energy can be written in the following 

way: 



F(T,P) 



dr 



^Tjexp 
r 2 



-n[T 2 +2M 



7=1 



01 



2ia'r 



n *3 

i=p—d+l 



01 



2i Rh 



X 



03(O\iT) 

e[(o\iT) 
e 2 (o\jT) 

6'Mir) 



0^ 



<9 4 (0| 



(3.11) 



where Vd is the volume in the non-compact directions parallel to the D p — D p system. 
This expression for the open string free energy will be our starting point in order to 
compute the phase transitions in the model under consideration. 



4 Phase transitions at finite temperature 

Given the explicit form of the effective potential eq. (13.2)) . it is interesting to see whether 
temperature corrections could modify the tachyon potential. We expect that at high 
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temperature the system is in a local minimum of the temperature-dependent part of eq. 
(13.21) . Then, as the temperature decreases, a point will be reached at which a second order 
phase transition will occur. The critical temperature 7^ for this to happen, as well as the 
relevant field space position T c can be found by solving the following set of equations: 

VZ ff (T e ,T e ) = and V? ff (T C7 T c ) = (4.1) 

where V e ff is given in eq. ( \3.2L and the / denotes d/dT c . 

In particular, in the case A = we expect that temperature corrections should lead to an 
effective potential in which the location of the minimum has shifted away from infinity. 
The physical reason for this is that moving towards T = can be thermodynamically fa- 
vorable: it costs energy, but it also reduces the mass of the tachyon and therefore increases 
the entropy of the tachyon gas [6], [7]. We will show explicitly that for temperature near 
the Hagedorn temperature the minimum will be shifted all the way to T = 0, in which 
case the open string vacuum would be stable. 

In the presence of a separation, it is interesting to see whether as the temperature de- 
creases the system will start rolling towards one or other of the zero temperature minima. 

4.1 Low Temperature 

As a warm up calculation and in order to check that our expression for the free energy of 
open strings, eq. (13.111) . reproduces known results in the limit of small separation between 
the Dp-i — D p _x pair (equivalently small A in the D p — D p T-dual system), let us study 
the low temperature approximation of eq. (13. lip . 

In |5], [7J it is shown that starting with the minimum of the potential (12. 8p at T = oo and 
at zero temperature, as the temperature increases the vacuum is shifted from T = oo to 
T = 0. In particular, it is shown that the position of the tachyon minimum, T m i n , moves 
almost linearly towards T = as the temperature increases. 

In this subsection we will recover this result in the more general background in which 
a Wilson line is present. 

In the large (3 limit, we can approximate the free energy (13.111) by the large r contributions 
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to the integral. In this limit the ^-functions become 

0i(O|zY) « 2e"^ 
6 2 {0\iT) « 2e"^ 
# 3 (0|zr) « 1 + 26"^ 

# 4 (0|zr) « l-2e^ T (4.2 
Using the above expressions, the free energy becomes 

-7T (T + 2a'A 2 (T - 4T 2 + 1) - 1) t - 7T f 



(4.3) 



F(T,/3) « -— ^ / drr"— exp 

This integral can be rewritten in terms of the modified K-Bessel function as 
where we defined 

f(T, A) = — + a 1 A 2 (T 4 - 4T 2 + l) 
In the limit in which both T and j3 are very large the free energy becomes 

F(T,(3) « -lZ¥L{2f{T,A))iexp (-M^/ff^j) (4.5) 

Inserting this expression in the effective potential (13.21) and minimizing it wrt T leads to 
the following condition 

TL i -^VmT~A) = (4.6) 

Ph 

In the case A = we simply have f(T, A = 0) = -| 2 and therefore we get: 

IU = ^ (4.7) 

If A ^ 0, let's assume that its absolute value is A < we are in the regime in which 

' V2ct 

T = is a maximum of the potential and the tachyon has negative mass near the origin. 
Then, if T is large but {A 2 T 2 ) < 1, we have 

Ph v 

and by expanding the square root 

T 2 mm - t 1 + a ' A2T ^ + T ^ w ^ 4 - 8 ) 

13 



The solution of the previous equation is either T m i n = 0, which is not in the T 3> 1 
approximation, or 

± ' mfal 4aM%/3 1 j 

Again here we assume that f3 is large, but (3 A <C 1 thus the square root can be expanded 
in terms of {(3 2 A 2 ) 



i 



/3„ 2*13 MVg „,,i> . 

r + = " — ^T" + ° ( ^ f ( 0) 

T_ is the solution we announced at the beginning of this section and it is in agreement 
with [6] and [7j. We see that, as the temperature increases, this minimum shifts almost 
linearly towards T — 0. |f| 

What is the meaning of the other solution, namely T + ? It has an opposite behaviour 
compared to the previous results: namely the minimum increases as the temperature in- 
creases. In fact, we see that this solution violates the approximation we made, namely 
A 2 T 2 « 1 and j3 A < 1. For example, say that A = 10" 7 , (3 « O(10 3 ) then we get 
T+ « 10 11 which gives A 2 T 2 » 1. 

At low temperature, no phase transition occurs regardless of the value of A. To see 
this, expand eq. (14.61) around large T, this time taking A = 0(1). Then 



dT V 



A (/3 H - 2y/2APiry/a F ) 



Computing the second derivative of the potential in this point, we have: 

«V = D^ = - ^ (4,2) 



OT 2 2-j23*A 

which is negative and clearly indicates the absence of a second order phase transition at 
low temperature. 

Finally, before moving on to consider the high temperature regime, notice that these 
results can also be found by considering the tachyon field alone, ignoring the contribution 



3 Note that the coefficient of the linear term in (3 in cq. (|4.9[) differs from [7] due to the different 
normalization we adopted in eq. p. lip . 



14 



of all other open string modes to the free energy of the system [6]. The reason is that as 
long as the temperature is low compared to the Hagedorn temperature, the tachyon has 
the lowest mass and its contribution is dominant. 

In this setup, the effective potential, e.g. of a Dg — D 9 system is given by the sum of the 
zero temperature tachyon potential and the free energy of the brane-antibrane system at 
finite temperature T = (3~ l . The one loop free energy density for the tachyonic degree of 
freedom in 9 + 1-dimensional space is given by 

T[Tj)= \S-Bf xo ^ i - p ~^ <4i3) 

where Uk = \j k 2 + M^ s , and Mns is given by eq. (I3.8P in which the bosonic degrees 
of freedom Nb and Nns are set to zero. Expanding the logarithm and performing the 
integration one gets: 

oo 

JF(T, /3) = - Y,(M^- 5 2-*M 5 NS K 5 (nPM NS ) (4.14) 

n=l 

where K 5 (z) is the modified Bessel function. At low temperature (large f3) we keep only 
the first term in the previous sum, obtaining 

F(T,f3) » -2(27ry 5 (M N s/P) 5 K 5 ((3M NS ) (4.15) 

which agrees with the free energy eq. (14. 5[) after we expand it in the limit in which both 
T and (3 are very large. 

4.2 High Temperature 

We will now use the expression for the free energy eq. ( 13.111) in order to investigate the 
behavior of the model at high temperature, that is to say, at a temperature close to, 
but below, the Hagedorn temperature. As we discussed earlier in section [3j in this case 
the canonical ensemble is generally not reliable and we should adopt the microcanonical 
description in order to compute thermodynamical quantities. But, as we argued there, for 
the case of D 9 — D 9 pairs with constant Wilson line, the canonical ensemble agrees with 
the computations made in the microcanonical ensemble. Thus we will focus our attention 
on the D 9 — Dg system with i^O. 

In contrast to the low temperature case discussed before, we now want to expand the 
integral in eq. (13.111) near r = 0. To facilitate this, it is convenient to introduce the 
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variable 



1 

r 



and consider the large t region expansion. 

Using the modular transformation of 9 functions and extracting the leading term in the 
large t region near the Hagedorn singularity, we obtain from (13.111) 



F(T,x) 



a 



2^ 2 /?r vn 



D 



i=p—d+l 



D+d-9 

dtt 2 exp 



p—d 

i=i 



iRjt 



n (T 2 + 2M 2 ) 



n * 

i=p— d+1 



— 7T (x 2 

ia't 



l)t 



x 



(4.16) 



-i^. In (14.161) we have in mind the case 



where A is a cutoff and we have defined x 
p = 9, d = 8, D = 1. 

We are interested in the behaviour of the system at the origin of field space, namely 
near T = 0, therefore, we expand the previous expression around this limit and we 
keep only the lower order terms. As shown in [8], the additional contributions from 
the quantized winding and momenta in (14.161) may modify the leading order Hagedorn 
singularity if the compactification radii are much bigger than the string scale. In the 
case where p = 9, d = 8, D = 1 there is only quantized momenta on the circle (since it 
necessarily lies in a direction parallel to the D9 ). If the radius of this circle is close to the 
string scale then as shown in [8] the Hagedorn singularity is dominant and the expression 
(I4.16P becomes 



F(T,x) 



cv p r 

■— J - / Gtrexp 
Ph J a 

p-9 



-ix (x 2 - l) t 



2a'A 2 7i 
t 



x 



f-¥ + (t 2 + 2a ' A 2 (T 4 - AT 2 )) 



where we have replaced p = d + D and defined 

C 



p—d 

V d \{Ri 



i=i 



(4.17) 



(4.18) 
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On the other hand the same will be true if we take the radius < vo? and assume that 
the energy of our system is sufficient to excite the quantized momentum modes along the 
S 1 . For the case p — 9,d = 8, D = 1 the latter condition means we may consider small Ri 
and large t such that R\t is still sufficiently large to allow us to approximate 6* 3 (o|^rj 
by unity 

Under this assumption, we can expand the exponential containing A in the previous 
expression, as long as A « 0(1). Keeping only the first two terms we find: 



F (T,x)^-^ [°° dte-^ 2 - 1 ^ \t E ^-7r(2a'A 2 (T A -4:T 2 + l)+T 2 )t E ^ 1 } (4.19) 

PH J A L - 1 

In the case in which p = 9 this integral can be easily done and the result is 



CV 

Ph 



7i (2a' A 2 (T 4 - 4T 2 + l) + T 2 ) T (0, tt (a: 2 - l) A) 



7T (X 2 — 1) 

(4.20) 

We may fix the cutoff scale A by comparison with the free energy computed in the micro- 
canonical ensemble with A = 0, [Tj. In particular, if we set A = (2tt)~ 1 , the two results 
agree. We have now all the ingredients to write the effective potential eq. (13.21) for the 
D 9 — Dg pair with constant Wilson line A, in order to study the phase transitions. By 
T-duality this is mapped to a separated D 8 — D 8 pair with separation d ~ \A\. 

The critical temperature (3 cr and the value of the tachyon T cr at which the phase 
transition occurs can be found by finding the solutions of the equation: 

= V = (4.21) 

ojv^m =v „ = o (4 22) 

We have: 

V' = T 9 e- T2 T(8A 2 a' -A(2A 2 a'T 2 + l)) + 

JL (2 (4T 3 - 8T) a' A 2 + 2T) V (0, n (x 2 - l) A) (4.23) 

PH 

A clear critical point is T cr = 0. Substituting this value in V"ff = gives the following 
condition 
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This equation is important, because it allows us to compute an approximated expression 
for the critical temperature at the point T cr = 0. However, we note that whereas the 
lhs of eq. (14.241) is positive definite, the rhs is positive definite only when < A < 4^ 
and A > A cr . Therefore there is no phase transition at T = for 4p < A < A c 



j. cr . 



When < A < 4f or A > A cr we can expand the gamma function in eq. fl4.24p near 
x = 1 using the fact that 

r(o,t) = - 7 -iogt + o(t) 

and we find 



Per ~ Ph 



16 (2A 2 a' 
1 + EXP <-^(8A^I)- 7 



(4.25) 



where we have set A = 1/27T. For weak coupling, g s needs to be small, therefore, the 
argument of the exponential in the previous equation is large and negative which means 
that the critical temperature is very close the the Hagedorn temperature. ( The limit of 
A — ► of this expression gives the results that Hotta found in [7].) 

Let us try now to find other solutions for the system of equations (14.211) and (I4.22p . 
Isolating the gamma-function from the first equation and substituting it into the second 
one gives the following condition for the presence of critical points: 

8e~ T2 T 2 T 9 (8 (T 4 — 3T 2 + 3) a' 2 A 4 + 2 (3T 2 — 4) a' A 2 + 1) 
A{T 2 -2)a'A 2 + l = ° (426) 

Except for the point T cr = 0, @ other possible solutions to the previous equation are 



2 _ 12a' 2 A 4 - 3a' A 2 ± a' A 2 yJ-ASa' 2 A i - 8a' A 2 + 1 

± ~ SAW ( 7) 

The argument of the square root is positive definite only for < A < 2 ^ 3a , , but for these 

values of A one can verify that T± is negative, giving an imaginary T. 

We conclude then that there is only a second order phase transition at T = 0. 



4.3 Phase Structure 

In order to study and understand the phase transitions within a system consisting of a 
Dg — Dg-pair at high temperature, standard thermal field theory reasoning can be very 

4 The divergence in the rhs of cq. (14.241) coming from the vanishing of the denominator is only apparent 
since it is due to the truncation to the second order term in the expansion of eq. (|4. 17)) around large t. 
The full expression of the free energy (|4.17p is not divergent for any value of A. 

5 The point T = +oo solves the equation (|4.26|) but it is out the range of our approximation, namely, 
we have expanded the free energy around T = 0. 
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useful: the minima of the effective potential at high temperature are located around those 
values of T which minimize the tachyon mass at zero temperature and hence increase the 
entropy of the tachyon gas. 

Recall that the mass of the tachyon in the presence of a Wilson line was given in eq. 
(I2.23P which we rewrite here for our convenience in terms of real T: 



M2 = ajm = L 

dd) 2 a' 



T — 1 



+ a' A 2 (T 4 - 4T 2 + l) 



(4.28) 



The extrema of M 2 (T) above are given by 

Ti = (4.29) 

and 



V8a' A 2 - 1 



(4.30) 



2AVa' 

The second derivative of eq. (I4.28P evaluated at T\ is (1 — 8a' A 2 ) /a' which is positive 
for A < \A cr . Therefore, for A < \A cr we expect that at high temperatures, T = is 
a minimum of the effective potential. If instead A > \A cri the point T = is a local 
maximum, the minimum being T2 ^ 0. 

We will now investigate the phase structure of our system in the 3 cases where 
< A < \A CT , 4f < A < A cr or A > A cr respectively. 

1 , 

4.3.1 < A < -A cr 

In this case, we know that there is a phase transition at T = which is also a minimum 
at high temperature. 

Referring to Figure^ we find that: 
1. When the temperature is slightly above to the critical temperature and close to the 
Hagedorn temperature, i.e. x ~ 1, we expect that the system to be in the minimum at 
T = 0. Therefore, for these temperatures and for < A < \A cr the open string vacuum 
is stable. 

We also see in Figure [2] that T = is actually a global minimum of the effective potential 
and the latter is negative around this point. To understand why this is the case, consider 
for example when A = 0. The zero temperature potential becomes V 1 {t=o,a=o} — 2T 9 
and the finite temperature contribution can be obtained from eq. f)4.20p 
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Figure 2: Phase Transition for temperatures close to the Hagedorn temperature and 
separation, < A < A cr /2. The plots show the effective potential as derived from eq. 
(I2.16P and eq. (14.161) using numerical integration, for various values of the temperature. 
We chose the values A = 0.3 A cr , g s = 0.1, a' = 1/2, A = 1/2tt, p = 9, d = 8, D = 1 in 
these plots. 

At the critical temperature, this expression can be rewritten using eq. (I4.25[) as 

F(0, x cr ) « -^exp fi?- - 7 ") (4.32) 

At weak coupling, g s is small and consequently the value of the free energy is much larger 
than the zero temperature piece, resulting in the effective potential becoming negative 
around T = 0. 

2. When the temperature is equal to the critical temperature given by (I4.24p the 
minimum at T = becomes flat and is uplifted so that the potential energy becomes 
positive. 

3. For temperatures lower than this critical temperature the point T = is a global 
maximum and the tachyon field will start rolling towards T = oo and the system will 
undergo tachyon condensation. 
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Moreover, we find that the value of the critical temperature for a second order phase 
transition at the point T = is proportional to the value of the Wilson line A: the greater 
the value of A, the closer the critical temperature j3 cr is to the Hagedorn temperature. It 
therefore requires more energy to produce a separated D§ — D 8 than a coincident one. 

4.3.2 4^ < A < A cr 

In this case, we know from eq. (I4.24p that there is no phase transition at T = 0, in fact 
there is no second order phase transition at all. At high temperatures, thermodynamic 
reasoning tells us that the system is in a global minimum of the effective potential which 
is located at T* « T2 given in eq. (14.301) while the point T = is a local maximum. Then, 
when the temperature decreases this minimum is uplifted and becomes shallower and 
shallower until it disappears at lower temperatures giving the zero temperature potential 
as the effective potential. (See Figure El) 




Figure 3: Phase Transition for temperatures close to the Hagedorn temperature and 
separation, 4^ < A < A cr . The plots show the effective potential as derived from eq. 
(I2.16P and eq. (I4.16P using numerical integration, for various values of the temperature 
and for the choice A = 0.7 A cr , g s = 0.1, a' = 1/2, A = l/2vr, p = 9, d = 8, D = 1. The 
temperatures are the following: in the left hand plot, (3 = 1.01 /3jj] in the right hand plot, 
(3 = (1.2, 1.3, 1.6, 2.8) (3h for the sequence of curves displayed from left to right. 



4.3.3 A > A cr 

In this case, referring to Figure |4] we find that: 

1. At temperatures very close to the Hagedorn temperature, the system is in a global and 
deep minimum of the effective potential, say T*, which is not at the origin of the tachyon 
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A > A cr A > A, 




Figure 4: Phase Transition for temperatures close to the Hagedorn temperature and 
separation, A > A cr . The plots show the effective potential as derived from eq. (12.161) and 
eq. (I4.16P using numerical integration and the following parameters choices : A = 1.1 A cr , 
g s = 0.1, a' = 1/2, A = 1/27T, p — 9, d — 8, D — 1. The temperature is the same in both 
plots: (3 = 1.01 (3h- The left plot is a zoom on the region of the effective potential close 
to T = whereas the right one shows the deep minimum at T = T*. 

space, T = 0. The point T = is a local maximum at this temperature, (see Figure HI) 

2. When the temperature approaches the critical temperature, given by eq. (I4.24p . 
the point T = becomes flat and there is a second order phase transition in this point, 
(see left plot of Figure [5]). However, the second minimum continues to exist and it is still 
deep. 

3. For temperatures below this critical temperature, we have two minima T = and 
T = T*. (See right plot of Figure \5§. As the temperature continues to decrease the 
minimum T* becomes shallow and eventually will disappear. 

4. Eventually, close to zero temperature, the minimum T* has disappeared and the 
system will undergo tachyon condensation. 

From our prospective, it seems unlikely that at zero temperature the system will be in 
the open string minimum T = but rather in the closed string minimum at T = oo. This is 
because, unless finite temperature tunneling effects happen between the two minima when 
the separation barrier is short, the system will likely find itself in the minimum T = T* at 
high temperature and as the temperature decreases, this minimum will become shallower 
and the tachyon field will eventually undergo tachyon condensation in the closed string 
vacuum. 
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Figure 5: Phase Transition for temperatures close to the Hagedorn temperature and 
separation, A > A cr . The plots show the effective potential as derived from eq. (12.160 and 
eq. (I4.16P using numerical integration, with the following parameter choices: A = lAA cr , 
g s = 0.1, a' = 1/2, A = 1/27T, p — 9, d = 8, D = 1. The temperatures are the followings: 
in the left plots (3 = (1.199, 2.4, )(3h for the lower and upper curves; in the right hand 
plot j3 = (7.8, 8.0, 8.6) f3jj for the three curves starting from the lower. 

5 Conclusions 

In this paper we have investigated the phase structure of a D p — D p pair at finite temper- 
ature, including a constant Wilson line A wrapping a spatial circle S 1 . By T-dualizing 
along the S 1 , this system is mapped to a D p _! — D p _! pair where the branes are parallel 
but separated by a distance d along the dual circle S 1 with d ~ \A\. Due to the limitations 
of the canonical ensemble as we take the temperatures close to the Hagerdorn transition, 
our results are mainly focused on the p = 9 case. The extension to all other values of p 
can be found by extending the microcanonical ensemble calculations of [7J, [8] with the 
inclusion of a non- vanishing Wilson line A. 

We found that the inclusion of A makes the effective potential acquire two minima at 
finite temperature if A > A cr compared to the situation with A < A cr (which includes 
the case A = of coincident D p — D p branes studied in [6] , [7] ) . This raised the question 
concerning which of the two minima our system is likely to be found. If we consider the 
case where we are at high temperatures, close to but below the Hagedorn temperature, 
then there is a single minimum with T^O, indicating the open string vacuum is unstable. 
As the temperature drops a second order phase transition occurs at the origin T = where 
a new minima develops which one can interpret as a meta-stable open string vacuum. 
However unless there are very special initial conditions it is unlikely that the system can 
be found in this metastable state but rather the second minimum at T ^ 0. The latter 
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coincides with the closed string vacuum T — > oo as the temperature approaches zero. 

If instead, A < A cr , then we showed that there is a phase transition in T = only for 
< A < 4p. In particular if this condition is satisfied, T = is a global minimum of the 
effective potential which is negative at high temperature and the system of a Dg — Dg-pair 
is stable. Then as the temperature decreases this minimum is uplifted and a second order 
phase transition occurs. 

Notice that in the dual picture, this implies that the separated D 8 — D 8 pairs undergo 
a phase transition, even in the case that the branes become coincident. This might appear 
at odds with the results of [7j, where it was found that no phase transition occurred in a 
coincident D p — D p pair with p < 9. But recall that in the dual system the D 8 — D 8 pair 
have one perpendicular spatial direction compactified on a circle. Thus the branes span 
all non-compact directions. In [8], phase transitions for a D p — D p pair were considered 
when some spatial dimensions are compactified on a torus. It was shown that a phase 
transition will occur for a coincident D p — D p pair even with p < 9 as long as the branes 
span all the non-compact directions. Thus our results are consistent with those in [8j. 

For 4^ < A < A cr we showed that there is no phase transition near T = 0. At 
high temperature the system is in a minimum away from the origin. As the temperature 
decreases this minimum eventually disappears. 

Our analysis in this paper is directed mainly at separated D 8 — D 8 pairs because of 
the limitations of the canonical ensemble for high temperatures. To properly study the 
case of separated D p — D p pairs with p < 7 (where we assume separation along a single 
compact direction) will require use of the microcanonical ensemble and extension of the 
complex temperature techniques used in the case of coincident D p — D p pairs, considered 
in [7J. Research in this direction is currently ongoing |21j . 

Finally, even if the D 8 — D 8 system were to find itself in the metastable minimum at 
T = as the temperature decreases, one should then consider the possibility that quan- 
tum tunneling effects can lead to the nucleation of closed string vacua T ^ 0. In the zero 
temperature case, [22] considered the possibility of tachyon tunnelling between the two 
minima of the effective potential when A > A cr . It would be interesting to extend this 
analysis in the case of finite temperature since then the barrier height and width between 
the two minima becomes a function of temperature so that it is not a priori obvious if 
tunnelling effects will be suppressed or not. 
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